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O ' Abstract 

o . 

CN . Cosmic strings form at the end of standard super symmetric hybrid inflation, 

> 

o 

' which is broken at the end of inflation is gauged B-L, there is a mixed scenario 

for leptogenesis: Right-handed neutrinos can be produced non-thermally during 



(N 

> 



reheating via inflaton decay as well as via cosmic string decay. We show that 
the parameter space consistent with CMB data can accommodate either or both 



(N . scenarios depending on the mass of the right-handed neutrinos. 

l> . 

^ : 1 Introduction 
o 

'-^ \ Supersymmetric theories beyond the standard model which predict light neutrino masses 

p • via the see-saw mechanism easily accommodate SUSY hybrid inflation PJ E] . Reheating 

^ I can proceed via inflaton decay into right-handed (s) neutrinos and thereby these models 

also provide an interesting framework for non-thermal leptogenesis jHHHE]. In the simplest 
• i-H , 

^ ■ version of hybrid inflation cosmic strings form at the end [ElIZ! ; if the inflaton sector couples 

■ 

I to right-handed (RH) (s)neutrinos these are B-L cosmic strings [HI EI- B-L cosmic strings 

are not superconducting (TUl E] . Most of the energy lost by the string network goes into 

gravitational radiation and right-handed neutrinos, and therefore these strings provide a 

second mechanism of non-thermal leptogenesis jH|. In this paper we calculate the relative 

contributions to the baryon asymmetry of the universe from reheating at the end of 

standard F-term inflation and from B-L cosmic string decay. 

Let Ggut denote a gauge group which contains the Standard Model gauge group 

as well as gauged B-L. F-term inflation requires the existence of a gauge singlet and 
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two Higgs superfields which transform in complex conjugate representations of Ggut! 
we assume that they break B-L when acquiring a non vanishing vacuum expectation 
value (VEV). In the case of SO(IO), they could be a 16 + 16 or 126 + T26. Note that 
the transformation properties of the Higgs representation can affect the stability of the 
strings; we shall not discuss it further here 0. Inflation takes place as the scalar singlet 
slowly rolls down a valley of local minima along which the VEV of the Higgs fields vanish. 
When the singlet falls below a certain critical value, the Higgs mass become tachyonic 
and inflation ends quickly in a phase transition during which the Higgs fields acquire a 
non-vanishing VEV breaking B-L spontaneously. If Ggut is semi-simple, the assumption 
of standard SUSY hybrid inflation then requires that Ggut breaks down to the Standard 
Model (SM) gauge group via at least one intermediate step, so that inflation solves the 
monopole problem; the gauge group broken at the end of inflation is not Gqut but an 
intermediate symmetry group Gint 3 U{1)b-l ^ U{1)b~l- During this phase 

transition, B-L cosmic strings form. The simplest example is Ggut = SO{10), E{6) or 
Pati-Salam, and Gint = SU{3)c x SU{2)l x U{1)r x U{\)b-l which breaks down to the 
SM gauge group at the end of inflation |T2l [H] • 

The GUT Higgs fields which trigger the end of inflation give heavy Majorana masses 
to the right-handed neutrinos. After inflation the universe is reheated by inflaton decay 
into RH neutrinos and sneutrinos. If the reheat temperature is less than the neutrino 
mass their out-of-equilibrium decay into (s)leptons and the SM Higgs(inos) produces a 
net lepton asymmetry [4^ . There is another contribution to the lepton asymmetry, coming 
from the decay of cosmic string loops 8J. Cosmic string loops decay into B-L Higgs and 
gauge fields, which in turn decay in right-handed neutrinos. Loops can also release right- 
handed neutrinos which are trapped as zero modes. 

Thus both inflation and cosmic strings contribute non-thermally to the baryon asym- 
metry of the universe. We investigate here which of these scenarios is most efficient using 
the observed CMB anisotropies as a constraint. The string tension, as well as the inflaton 
mass are set by the symmetry breaking scale at the end of inflation rj (the B-L breaking 
scale) and the Higgs self quartic coupling k (the superpotential coupling). The string 
tension also depends logarithmically on the gauge coupling constant which we set to the 
unification value for the MSSM. The resulting lepton asymmetry depends on r], k and 
the RH neutrino masses. Requiring that the inflaton gives the observed density pertur- 
bations fixes r/ as a function of n. The string contribution to the density perturbations is 
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constrained to be less than about 10% [Ij. This restricts the value of as discussed in 
our previous paper [T5] . 

Most of the lepton asymmetry in B-L string decay is generated at the earliest time 
at which leptogenesis is possible. Refs. [ini E] assume that the reheat temperature is 
high enough for the lightest RH neutrino with mass Mi to be initially in thermal equilib- 
rium. Any asymmetry generated is washed-out until the lightest RH neutrino freezes out. 
Hence, in their analysis the lepton asymmetry is dominated by the contribution generated 
at the freeze-out temperature T ~ Mi. However, there is also the possibility that the 
RH neutrinos are never in thermal equilibrium after inflation. Then the loops formed 
immediately at the end of inflation give the dominant contribution to the lepton asymme- 
try. The final asymmetry is determined by three factors. First, it depends on the initial 
string density, which can be different from the density during the scaling regime. Second, 
the universe goes from matter domination to radiation domination during reheating. The 
earlier this happens, the larger the asymmetry. And the third factor which plays a role is 
the CP asymmetry per decaying (s)neutrino, which depends on the details of the neutrino 
sector. 

The paper is organised as follows. In Sec. 2 we review standard SUSY F-infiation 
coupled to N=l SUGRA and discuss the CMB constraints ^H]- In Sec. 3 we determine 
the parameter space for successful non-thermal leptogenesis from reheating at the end of 
infiation. In Sec. 4 we turn to the parameter space for successful non-thermal leptogenesis 
which results from the decay of B-L cosmic strings. We distinguish three different cases. 
The first possibility is that Mi < and Mi > Tr, with Mi the lightest RH neutrino 

mass, the infiaton mass and Tr the reheat temperature. In this case reheating goes via 
production of RH neutrinos which are out-of-equilibrium at the end of inflation, and both 
non-thermal leptogenesis scenarios compete. The second possibility is that Mi > 
and reheating is gravitational. Then cosmic strings give the sole contribution to the lepton 
asymmetry. Finally there is the possibility that Mi < Tr and the RH neutrinos are in 
thermal equilibrium at production. In this case both cosmic string decay and standard 
thermal leptogenesis contribute to the lepton asymmetry. In Sec. 5 we consider the 
consequences on leptogenesis of generating a dynamical /x-term. We give our conclusions 
in Sec. 6. 
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2 Hybrid inflation & CMB constraints 



In this section we summarise the bounds on the parameter space imphed by current CMB 
data for standard hybrid inflation with cosmic strings. The details can be found in Ref. 

m 

The superpotential for standard hybrid inflation is E] 

Wm = kW-^'), (1) 

with S a gauge singlet superfield, and (j), cj) Higgs superfields in AT-dimensional complex 
conjugate representations of a gauge group Ggut- Upon acquiring a VEV the Higgs fields 
break Gint 3 U{1)b-l down to a subgroup H 7^ U{1)b-l- The supersymmetric part 
of the scalar potential is given by (we represent the scalar components with the same 
symbols as the superfields) 

^susY = - r/T + ^'\S\'m' + |<^n + Vd. (2) 

Vanishing of the D-terms enforces |0| = \(f)\. Assuming chaotic initial conditions the fields 
get trapped in the inflationary valley of local minima at \S\ > Sc = rj and = = 0. 
The potential is dominated by a constant term Vq = k'^i]^ which drives inflation. Inflation 
ends when the inflaton drops below its critical value Sc (or when the second slow-roll 
parameter rj equals unity, whatever happens first) and the fields roll toward the global 
SUSY minima of the potential |0| = \(f)\ = rj and S = 0. During this phase transition B-L 
cosmic strings form [Bl H] • For a discussion on various GUT models see Ref. [HI . 

The scalar potential in Eq. ^ gets corrections from SUSY breaking by the finite 
energy density in the universe during inflation (given by the Coleman- Weinberg for- 
mula), from SUSY breaking today, and from supergravity. The hidden sector expec- 
tation values responsible for low energy SUSY breaking can generically be written as 
(z) = arup, (VFhid) = fJ'fnl, (^^^7^) = cfirup, with z a hidden sector field, a, c dimen- 
sionless numbers, and fi a mass parameter related to the gravitino mass via 1713/2 = e'^l^/^/i. 
Further rUp = (SvrG)^^/^ = 2.4 x 10^*^ GeV is the reduced Planck mass. The scalar poten- 
tial along the inflationary valley can be calculated using the SUGRA formula 
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Assuming a minimal Kahler potential, the scalar potential including all corrections is 



2^4 



V = K, 1] 



+ 



K 



327r2 



2 In +iz + 1)2 ln(l + z^') + (z - if ln(l 



|a|V 
3mp 2mp 



+ Ky4m3/2r7 a, (4) 



with a = \S\/\^ the normalized real field, A = 2v^ cos(arg/i — argS*), and we assume 
that argS* is constant during inflation. The cosmological constant today vanishes for 
|c + a*p = 3, and we have dropped subdominant terms. Further we used the notation 
z = x"^ = IS'p/r/^ = a'^/{27f) so that z = x = 1 when a = a^. The flrst line is the tree 
level potential term, the second line is the one-loop Coleman- Weinberg correction due to 
SUSY breaking during inflation 0, and the third line are the SUGRA corrections. 

The Coleman- Weinberg potential and the non-renormalisable terms are always present, 
independent of low energy SUSY breaking. The A- and mass terms can be made small 
for a small gravitino mass (as in gauge mediation), or by tuning A and/or a. Another 
possibility is that the hidden sector superfleld z only acquires its VEV after inflation, so 
that these terms are absent during inflation. We note that, assuming gravity mediated 
SUSY breaking, the generic values ms/2 ~ lO^GeV, A ~ 1, a ~ 1, give rise to too large 
A- and mass terms, incompatible with the CMB data |15j. 

Both strings and the inflaton contribute to the primordial density perturbations [^ll^j. 
Cosmic strings do not predict the measured acoustic peaks in the CMB and hence their 
contribution to the temperature fluctuations should be small, less than about 10% T^. 
The string contribution is proportional to the string tension 

fx = 27rr]'e{P), (5) 

with f3 = (m^/m^)2. The Higgs mass is = k^i]'^ and the vector boson mass is 
m\ ~ QQUT^ff with the GUT coupling Qq^'^ ~ 47r/25. When /3 = 1, the strings satisfy 
the Bogomolny bound (this is the case of cosmic strings which form at the end of brane 
inflation) and ^^(1) = 1. However in the case of SUSY GUTs, the strings never satisfy the 
Bogomolny bound, [3 <1 always, and 



1.04/?o■l9^ /3 > 10-2 
(3 < 10" 



2.4 ,.,^_2 (6) 



log(2//5) ' 
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Requiring the string contribution to the quadrupole to be less than 10% gives the bound 

m 

Gii < 6.9 X 10-^ r/bnd < 4.1 x lO^^^j^M., (7) 

Here y parameterizes the density of the string network, and should be taken form numer- 
ical simulations. Recent work predicts y = 9 ±2.5 jTH]- Older simulations give y = 6 PU] . 
and semi-analytic approximations give ?/ = 3 — 6 j2I] . 

The density perturbations produced in hybrid inflation can be calculated using the 
slow roll formalism for the potential in Eq. (j3)). Setting it equal to the value observed by 
WMAP gives as a function of k. We will use the analytic approximations, derived in 
the limit that the Coleman- Weinberg(CW)-potential respectively the non-renormalisable 
(NR) terms dominate the potential: |15j 

rj,^ = 5xl0i^GeVArV3(_^y/\ (8) 
r/NR = 3 X 10^^ GeV (^) . (9) 

The symmetry breaking scale is restricted to the range 

Vcw ^V< min[?7NR, ^?bnd]- (10) 

If the A- and mass terms are absent or subdominant during inflation there are two distinct 
solutions, corresponding to r]^^ and t^nr (the upper bound T^bnd comes from the fact that 
the string contribution to the CMB is limited). If on the other hand these terms do play a 
role, the whole range is possible. This is illustrated in Fig. ^ which shows 77 as a function 
of K, for A/" = 1, 16, 126. In this plot the A- and mass term are assumed to be negligible; if 
they are not the solution is somewhere in the range given by Eq. The straight part 

at relatively large coupling is well approximated by r]^^. At low value there is a second 
branch of solutions given by ?7nr. Also plotted is r^bnd for y = 3; above this line the string 
contribution to the CMB is more than 10%. The CMB constraints are satisfied for the 
coupling range 10"^ ^ ^ lO^^/A/" and the SSB scale range rj ~ 10^^ — 10^^ GeV. 

The CMB bound can be avoided if the strings are semi-local or not topologically stable 
down to low energy and decay at some later phase transition [0]. 
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Figure 1: t] vs. k for A/" = 1, 16, 126 and y = 3. Further shown is the 10%-bound. 

3 Non-thermal leptogenesis from inflaton decay 

In this section we review the non-thermal leptogenesis scenario which happens during 
reheating as a result of inflaton decay into right-handed (s)neutrinos [HE]. 

The Higgs fields (p and (p break local B-L spontaneously when developing a VEV at 
the end of inflation. The right-handed neutrinos acquire super-heavy Majorana masses 
via their coupling to (p. Some other GUT superfield (this is model dependent) gives 
a Dirac mass to the neutrinos and the light neutrinos acquire a super light Majorana 
mass via the see-saw mechanism j22I- The lepton asymmetry, generated as the right- 
handed (s)neutrinos decay into SM Higgs(inos) and (s)leptons, is converted into a baryon 
asymmetry via sphaleron transitions 

Depending on the transformation properties of the Higgs representation, the right- 
handed neutrino masses are generated via normalisable or non-renormalisable superpo- 
tential terms 



W 



— 'yij(f)(f)FiFj, 



(11) 



W 



yijCpFiFj, 



(12) 
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where F is the n-dimensional spinorial representation of G which contains the right- 
handed neutrino superfield A^, and i,j = 1..3 for three famihes. At the end of inflation S 
and 0+ = (50 + 50)/v^, with (p = r] + 5(l) and (p = r] + S(j), oscillate around the global SUSY 
minimum of the potential until they decay into right-handed neutrinos and sneutrinos, 
thereby reheating the universe |1] . We work in the basis where the right-handed neutrinos 
mass matrix is diagonal. The decay rates r(0+ — >• iYjiVj) and T{S NiNi) are equal and 
given by 

2 

with X = S*, the oscillating fields which have equal mass = kt], and Mj = ?/j?7, 
the mass of the heaviest right-handed neutrino Ni (sneutrino Ni) the infiaton can decay 
into (i.e., which satisfies Mj < m^/2). The reheat temperature is then 

where we have used g^, = 228.75 for the MSSM spectrum. 

Non-thermal leptogenesis from infiaton decay takes place if the following constraints 
are satisfied: 

Kinematic constraint: Infiaton decay into right-handed (s)neutrinos is kinematically 
allowed if Mj < i.e., 

Mi{K) < -K7] (15) 

Gravitino constraint: Gravitino overproduction is avoided if the reheat temperature 
Tr < 10^° GeV i22j. This gives 

M.W<1.6xlO"Gev(^)^. (16) 

The upper bound on Tr is model dependent and can be as low as 10^ GeV. We 
note that the gravitino constraint can be avoided if the gravitino mass is sufficiently 
large so that it decays before BBN. 

Gravitational decay: The decay rate into right-handed (s)neutrinos should be larger 
than the gravitational decay rate into light particles. In a full theory, the superpo- 
tential is 1^ = VFinfi + Wud + W^gut, where W^gut contains GUT superfields, some 
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of them containing the MSSM fields. The gravitational decay rate of the infiaton 
into light SM particles can then be computed by considering for example a term of 
the form VF^gut = aHFF, F containing the standard model fermions and H some 
GUT Higgs superfield containing the SM Higgs. In the SUGRA potential Eq. Q 
there is a coupling between the infiaton and the SM particles, leading to a decay 
rate 

1 mr,?7 



X 

i7T 



fgrav — „ 4 • (17) 



Note that this is parametrically smaller than the standard gravitational decay rate 
Tgrav = (l/87r)m^/mp 24^. Requiring > Tgj-av then leads to 

^ #• (18) 

The wash-out constraint: The lepton asymmetry produced by the decay of the RH 
neutrino Ni is washed out by the L-violating processes involving RH neutrinos, 
unless they are out of thermal equilibrium which is automatic if Mj < Tr, with 
j = 1, 2, 3. The strongest constraint is for the lightest RH neutrino: 



^'^^-> > 1 ,/!^, (19) 



Mi{K) ~ ley ri{K) 

with as before the mass of the RH neutrino the infiaton decays into. This 
implies Mi > 10"^ - IQ-^Mi for (r^, k) allowed by CMB data (see Eqs. l8llTn|l. No 
wash-out is assured if Mi < < M2, M3 and the infiaton decays into the lightest 
right-handed neutrino. The CP violating parameter e can be improved by at most 
a factor (M2/M1) ~ 10^ — 10^ if the decay is not into the lightest neutrino, but into 
the next to lightest one (see the Appendix: ei oc Mi and 62 oc M2) [Sj. 

Perturbative couplings: We require the couplings 7^, yi in Eqs. fllH IT^ to be less 
than unity. For a renormalisable mass term this bound cannot compete with the 
kinematic constraint. For a non-renormalisable mass term this implies 

2 

Mi < ^ ~ 1 X 10^^ GeV(Arft)2/3, (20) 

where in the last step we used t] = rj^^, see Eq. ((HI). If rj/rrip < k, which only 
happens for k > 10~^, then for non-renormalisable mass terms all neutrino masses 
are lighter than the infiaton mass. 
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Lepton asymmetry: The lepton asymmetry produced is 

ut 3 Tt) , . 

— = 7; — ^^, 21 

s 2 

with the CP asymmetry per decaying RH neutrino Ni. For hierarchical RH 
neutrino masses and hierarchical light neutrinos the CP asymmetry in the decay of 
the lightest RH neutrino is bounded by [201 121] 

|,,| < 2 X 10-- (^) . (22) 

' ' - VlOeCeVy V 0.05 eV J ^ ' 

As discussed in the appendix, the upper bound on €2 is of the same order of mag- 
nitude, but with Ml — > M2 in the above formula. The CP-asymmetry induced 
by the decay of the heaviest RH neutrino is suppressed by a factor M2/M3. For 
quasi-degenerate light neutrinos (mi ^ m2 ~ ^ (Am^^-jj^)^/^) the asymmetry is 
smaller [21] 

,nf Mi \ /(Am2 )i/2\ /(Am2 )i/2\ 
with m = I/Sa/wiI + + m|. 

If the RH neutrinos are quasi-degenerate and Mi — Mj ~ Fj the CP-asymmetry is 
enhanced. The only constraint is then |ej| < 1. Upper bounds on the CP-asymmetry 
in type II see-saw models have also been derived for non-degenerate neutrinos, and 
are of the same magnitude as Eq. (j22|) |28] . 

The baryon asymmetry inferred from BBN translates into a primordial lepton asym- 
metry given by ni/s = 2.4 x 10~^° for the MSSM spectrum. For hierarchical light 
neutrinos, this is obtained for 

Mi{k) > 5.4 X 10=^ f^fli!^\ ' /GiV. (24) 
V "^p / 

For smaller Mj the produced asymmetry is too small. 

The parameter space compatible with NT leptogenesis from inflaton decay is shown in 
Figs. 121 and El The kinematic, perturbative coupling and gravitino constraints all give an 
upper bound on M,. The perturbative coupling constraint is weakest and not shown. The 
kinematic constraint is strongest for small ft, whereas the gravitino constraint dominates 
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Figure 2: Mj vs. k for r] = rj^^ and A/" = 1. The parameter space is bounded by the 
kinematic constraint (left), successful leptogenesis (right), and the gravitino constraint 
(lines, for Tr = 10^ 10^ 10^ 10^° GeV) . 



for large coupling. The gravitational decay and leptogenesis constraint give a lower bound 
on Mj. The leptogenesis bound is strongest. The k range is bounded by the CMB data, 
as given by Eqs. (0 El EII) and Fig. d to 10^*^ <k< 10~^/Af. 

Fig. 121 shows the parameter space for rj = t]^^ and JV = 1. The colored regions are 
excluded. The upper bound is fixed by the reheat temperature; the bounds for Tr = 
10^, 10^, 10^, 10^° GeV are shown. Leptogenesis is only possible for Tr > 10^ GeV. For 
Tr ~ 10^ GeV, successful leptogenesis requires Mj = 10^ — 10^^ GeV and k = 10~^ — 10~^. 
One should remember that in addition to the constraints shown in the plot, it should 
be checked that Mi > Tr and there is no wash-out of asymmetry. This is the case for 
Ml > 10~^Mj, with leptogenesis dominated by inflaton decay into Ni, in agreement with 
Eq. (HH). 

For a fixed coupling value, the possible range of neutrino mass Mj is about a decade. 
Therefore the bounds are all close to saturation. E.g. the leptogenesis constraint gives 
Mi (X (ei/ef"^"")"^/^, and thus (ej/ef"^'') > 10^^ is required. Hence, degenerate light neu- 
trinos with rh ~ eV are marginally excluded. 
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Figure 3: Mj vs. n for max[r7Q^, //nr] < r] < r/bnd, = lO^'^^GeV and Af = 1,16,126. 
The parameter space is bounded by the kinematic constraint (left), successful leptogenesis 
(bottom), the gravitino constraint (top), and for J\f = 16, 126 by the CMB data (right). 

Fig.Olshows the parameter space for t] in the whole range of Eq. (fTn|) for Tr = 10^° GeV 
and A/" = 1, 16, 126. The parameter space is enhanced compared to Fig. El about twice as 
big. The main consequence of increasing A/" is that the small k range is excluded by the 
CMB data, and that leptogenesis requires a slightly larger neutrino mass. 



4 Leptogenesis from string decay 

Cosmic strings form in SUSY GUT models with standard hybrid inflation jHl EI- The 
string mass per unit length is then constrained by CMB data, see Eqs.dS}, & and 
When the symmetry broken is gauged B-L, they also provide a non-thermal scenario for 
leptogenesis [H|. We first describe various possible NT leptogenesis scenarios with B-L 
strings forming at the end of inflation. We then discuss the evolution of the string network 
and analyse the various scenarios in details. 
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4.1 Leptogenesis scenarios 

Depending on whether the mass of the hghtest RH neutrino Mi is larger or smaller than 
the reheat temperature Tr, and on whether the inflaton decays into right handed neutrinos 
or not, there will also be a contribution to the lepton asymmetry of the universe from 
inflaton decay or standard thermal leptogenesis. We distinguish the following cases: 

Case 1: The reheat temperature is lower than the lightest right-handed neutrino mass, 
Tr < Ml, and there is no wash-out at any time. The lepton asymmetry is set 
by the earliest time that string loops form, which is right at the end of inflation. 
Reheating of the universe takes place at a later time, via inflaton decay into right- 
handed neutrinos. Apart from NT leptogenesis from B-L strings, there is also 
a contribution from non-thermal leptogenesis from inflaton decay as discussed in 
section El 

Constraints: gravitino constraint Tr < 10^'^ GeV, gravitational decay constraint 
Tgrav < '^Ni kinematical constraint Mj < 

Case 2: Same as case 1 but now the inflaton does not decay into RH neutrinos. For 
example, if Mj > my./ 2 \/i and there is no other superpotential term involving the 
singlet or the Higgs flelds, decay is through gravitational interactions: Fat < Fg^av 
Gravitational reheating can alleviate the gravitino constraint. Leptogenesis comes 
from the decay of the string forming gauge fleld in right handed neutrinos. This is 
the only contribution to the lepton asymmetry. 

Constraint: gravitino Tr < 10^'' GeV. 

Case 3: The reheat temperature is higher than the lightest neutrino mass Tr > Mi. 
The asymmetry will be washed out at high temperatures by L-violating processes 
mediated by A'^i, and can only be created for Tr < Mi. The asymmetry is then 
dominated by the loop formation rate at T ~ Mi. There are two contributions 
to the lepton asymmetry: NT leptogenesis from B-L strings and standard thermal 
leptogenesis. This is the case considered before jl6| I17j. 

Constraints: gravitino constraint Tr < 10^° GeV, kinematical constraint Mi < 

In order to analyse the various scenarios, we study analytically the evolution of the 
string network, both in the scaling regime and at the initial times. 
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4.2 String network and neutrino density 

The evolution of a cosmic string network has been extensively studied over the years 
|29j . Numerical simulations and analytical studies agree that the string network reaches 
a scaling regime, in which the energy-density carried by the network remains a constant 
fraction of the total energy density in the universe. The scaling solution is an attractor 
solution, and is independent of the initial string density. This is one of the reasons that 
string network at formation has not been discussed much in the literature. However, the 
lepton asymmetry is dominated by the initial time, and thus depends sensitively on the 
initial density. We first discuss the familiar scaling regime, before turning to a discussion 
of the initial string density. 



4.2.1 The scaling regime 

To describe the approach to the scaling regime, we introduce the characteristic length-scale 
L which sets the correlation length and the average distance between long strings [5^. 
The energy density in long strings, poo ~ /w/-^^, evolves as 

p^ = -2Hp^-f{p)^ (25) 

where the terms on the right hand side describe the energy loss due to expansion of 
universe, and due to production of loops respectively. The function f{p) depends on the 
reconnection probability p as f{p) ~ ^/p. For gauge field theory cosmic strings p = 1. 
Introducing 7(t) such that 

L = ^{t)t, (26) 

one finds that the above equation has a stable attractor solution, the scaling solution. It 
does not depend on the initial string density: 

where we wrote P = Ht. Since 7^ = y/p/{2{l — [3)) = 0{1) is constant, the long strings 
scale with the horizon. 

The scale L characterizes the network on macroscopic scales, but does not say anything 
about what happens on the smallest scales. Simulations show that the long strings have 
small-scale wiggles, whose characteristic length also scales with time These wiggles 
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set the typical loop size, which we parametrise as 



hoop ~ at (28) 



with a ~ (FG/i)", and T ~ 50. The "standard" value is n = 1 giving a = ai = (FG/i) (4 
More recent simulations suggest n = 3/2(5/2) in the radiation (matter) dominated era 
j3Uj . The loop formation rate nioop is set by the requirement that the string network 
keeps its scaling solution. The loops loose energy by emitting gravitational radiation and 
contract until the loop radius becomes of the order of the string width, at which point 
it decays emitting X-particles (with X = (f),A,N, i.e., the string Higgs or gauge fields, 
or RH neutrino zero modes). For a < ai, the loop lifetime is less than a Hubble time 
and we can then neglect the red shifting between birth and death. The injection rate of 
right-handed neutrinos during scaling is then simply 

riN = XNTlloop - ^^f^, (29) 

where xn is the number of right-handed (s)neutrinos produced per decaying loop. 

The minimal number of RH neutrinos released per loop is xn = 1- However, we expect 
the loop to decay when its radius becomes of the order of the string width m^^ with a 
burst of the strings Higgs and gauge particles [23 ■ These in turn (mostly) decay into RH 
neutrinos. The number of neutrinos emitted per loop is then of the order 

with ^(/?) ~ 0.1 — 1 given in Eq. (jH)). If the loops decay mainly into Higgs fields mx = fri^ 
and xat ~ k^^, whereas if decay is mainly into gauge bosons mx = m^ and x^ ~ n^^. 
If the string width at which the loop decays is smaller than the the inverse Higgs mass 
m^^, then a; a? is correspondingly smaller. 

The lepton density is obtained by integrating Eq. with a red shift factor (a(tin) / ci{t))^ 
{tin/tf'^'^ to account for the expansion of the universe. Here it is assumed that both t and 
the initial time tin are in the radiation dominated era following inflaton decay. This gives 
(using t ~ H~^) 

n,{H) ^ ^i/f (31) 
independent of the reheating temperature. 
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Some simulations find that loops form on the smallest possible scale (given by the 
resolution of the simulation). These results suggest that scaling is maintained mostly by 
particle emission rather than via loop formation and subsequent gravitational decay. This 
is called the VHS scenario, after the authors of Ref. [211 . The emitted X-particles are 
the RH neutrinos themselves, due to the existence of zero mode solutions, and the string 
forming Higgs and gauge fields which then (mostly) decay into RH neutrinos. This gives 
a RH neutrino injection rate jHSl 

hN^f^-^,t-'. (32) 

with X = (f),A, N, and the fraction of the energy in loops that goes into X-particles. 
Not all of the loop energy can go into high-energy particles, as this would give too large a 
diffuse 7-ray back ground, in conflict with the EGRET data |3S] ■ Combining the EGRET 
bound with the CMB bound gives^ 



We note, however, that the EGRET flux is dominated by late times. It is therefore not 
impossible that is time-dependent, and much larger than the bound above at early 
times. 

When leptogenesis takes place after reheating, the lepton number density is obtained 
by integrating ()H2j) taking into account the expansion of the universe. This gives 

riLiH) ^ (34) 

7s 

4.2.2 The initial string density 

Cosmic strings are formed during the B-L breaking phase transition. The string density 
is set by the correlation length at the time of the phase transition ^ |H4| IHH] . The universe 
is cold at the end of inflation, and the equilibrium correlation length is set by the mass of 
the symmetry breaking Higgs fields ^(t)~^ = rn^{t) = K,'^{S'^{t) —rf). Both the correlation 
length and the relaxation time r = ^ diverge during the phase transition, and eventually 



similar injection rate and EGRET constraint apply for large loops which undergo 'quick death', 
i.e., loops that decay through many self-intersections into small loops which in turn decay emitting heavy 
X-particles. 
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4> must fall out of equilibrium. The correlation length at freeze out ^ is thus in the range 

(nr^)-' < i < H;\ (35) 

with H^: the Hubble constant at the end of inflation. The lower bound is set by the 
maximum Higgs mass = krj obtained in the vacuum. The upper bound is set by 
causality, as fluctuations cannot exceed the horizon. A more careful estimate of ^ is given 
in Ref. [^E!- Writing S(t) = Sc — St near the phase transition, the inverse Higgs mass, 
which sets the correlation length, is Tn'^{t) = —{K'^riS)t. Freeze out happens approximately 
at the time when the relaxation time is equal to and thus 

i ^ i{-r) = {k'vS)-'/\ (36) 

During slow roll inflation 5* ~ {60H^). If the velocity does not change much between 
time observable scales leave the horizon and the time of the phase transition, then this 
is a good estimate. The freeze-out correlation length ^ determines the typical distance 
between cosmic strings. 

At the time of string formation (quantities will be denoted by subscript the corre- 
lation length ^ = L = •j^^t^, (see Eq. fj26|l ) 

r^H^r^^. (37) 

TUp 

Using Eqs. (fH^ IH^ . we get the range and "best" value for 7*: 

7mm = ^?Mp < 7* < 1 = 7max, 7bcst = 7* ~ 0. (?7/mp) (38) 

Since 7^, < 7^ and poo oc 7~^, the energy density in the network is initially larger than 
during the scaling regime. Solving Eq. (j25|l with the initial conditions above, it can be 
seen that the scaling regime is typically reached in only a couple of Hubble times. We will 
use the approximation that initial network reaches the scaling regime instantaneously, 
and an amount of energy ~ ji/{j^t^Y is dumped into loops and/or particles at the 
initial time given by Eq. ()37|1 . ^ 

Eq. (j^H|l for the loop size breaks down at the initial time. The reason is that the loop 
radius is smaller than ~ ""^I^, ('^a^)' width of the profile function of the Higgs 



^We neglect friction in the thermal bath which is absent for T < G^m-p. Note that the strings 
form before reheating has completed. Taking for T the reheat temperature Eq. (|14|) . this gives Mi < 
^Qjfl'^ I^K^^'^mJ'^) - 6 X Ge\eK^/^U^/^, where in the last step we used r/ = Ty^w 
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(gauge) field, and various parts of tlie loop overlap. It is not possible to speak of cosmic 
string loops anymore, which are well defined only for /loop ~ at > 27rm^^, (27rm^^), which 
requires t ^ t*. Therefore it is expected that initially the energy loss in gravitational 
radiation is small, and the network mainly decays directly into RH neutrinos and into 
gauge and Higgs fields which (mostly) decay into RH neutrinos. The lepton number 
density ai H < Tr is then 



per decaying RH neutrino, see Eq. ()22|) . Before inflaton decay the universe is matter 
dominated and a oc t^/^, afterwards the universe is radiation dominated and a oc t^^"^. 
This is used in the second step to write (a(t*)/a(t))^ = H^^'^H^'^T^^'^ , where the time 
t ~ is after reheating of the universe H < T. This factor takes the red shift due 
to the expansion of the universe into account. The earlier the transition from matter to 
radiation domination, i.e., the larger the decay rate, the larger Tr and the larger is the 
final number density. The initial time tin is right at the end of inflation, see Eq. (jHTjl . 

Eq. (j39|) is the same as the lepton asymmetry produced during the scaling regime as 
given by Eq. (jH^ under the replacement 7* — > 7s (difference in energy densities stored in 
long strings), and F (difference in whether leptogenesis takes place before or after 

reheating) . 

4.3 Results 

For the parameters at hand T < and the inflaton decays some time after inflation. 
The lepton asymmetry ni/s is to be evaluated after reheating, when entropy is defined. 
The entropy is s = {2tt^ /A5)g^T^ ~ lO^T^ with T ^ OA^/rri~H. 

4.3.1 Case 1 

The string contribution is dominated by the asymmetry produced at the initial time. The 
lepton number density is given by Eq. p9|l . Dividing by the entropy gives 




(39) 



where is the fraction of the energy going into X-particles and is the CP asymmetry 



■5 7in mxm\ 



3/2 ■ 



(40) 
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Using Eqs. (fT^ and (j221)) we see that n^/s is proportional to MiMj the mass of the 
heaviest RH neutrino the inflaton can decay into and the mass of the RH neutrino which 
is mostly produced by strings. For Mj = Mj, which is automatic if the strings decay 
mostly into 0-particles, successful leptogenesis with rii/s = 2.4 x 10^^° requires 

where in the second step we used r] = r]^^, and we introduced 



C is minimized for 7^, 'jmin and — ^ 1. In this limit the energy density in the string 
network is of the same order as the energy density in the oscillating inflaton field. If both 
the inflaton and the strings decay into (p particles {X = (p) this gives a contribution to the 
lepton asymmetry of similar magnitude. If however 7* > 7min or < 1, the energy stored 
in the string network is subdominant, and thus also its contribution to the asymmetry: 

('^L)infiaton f ( ^* ] ('^3) 



strings 

The string contribution can dominate over the inflaton contribution to the asym- 
metry if the strings decay mostly into RH (s) neutrinos. Another possibility is that 
Ml < < M2 and the string decays mostly into gauge particles. The inflaton de- 
cays into lightest RH (s)neutrino, whereas the gauge field can also decay in the next to 
lightest one. Since £2/61 ~ M2/M1 the CP-asymmetry is then larger per decaying A- 
particle than per decaying 0-particle. However, the number of gauge particles produced 
by string decay can be of the same order as the number of Higgs particles from inflaton 
decay only in the limit n —>■ 1 (so that ~ 171^) and — > 1. 

For degenerate light neutrinos the CP asymmetry is smaller by a factor Amatm/'"^3, 

1 /9 

and Mi oc e/ is larger. The A/"- dependence of the neutrino mass Mj is weak; the main 
effect of considering a larger Higgs representation is the stronger constraint on the coupling 
coming from CMB data: k, < lO^^/A/". The bound on the neutrino mass is weakest in the 
limit K, 1. This is only possible if the strings do not contribute to the CMB anisotropics, 
i.e., if they are semi- local or are not topologically stable jHj. 
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Figure 4: Mj vs. k for case 2 with rj = t]^^ and A/" = 1. The hnes A-C correspond 
to A= (7in,X, = (7min,^, 1), 5 = (7best,^, 1), C = (7min, v4, 10"^) . The parameter 
space is bounded by the kinematic constraint (top and bottom) , and gravitino constraint 
(parallel lines, for Tr = 10\ 10^^ GeV). 



4.3.2 Case 2 

If < Mj < 7] inflaton decay into RH neutrinos is kinematically forbidden, and in 
the absence of other direct couplings the inflaton decays gravitationally. This implies a 
low reheat temperature which alleviates the gravitino constraint. Note however that the 
lepton asymmetry ul (x F^/^ is less efficient and large neutrino masses are needed. The 
string loops decay into RH neutrinos and gauge quanta which subsequently decay into 
RH neutrinos. 

To get the RH neutrino mass required for leptogenesis in case 2, we use Eq. (|in|) with 
Ftv replaced by the gravitational decay rate, to yield 



M - 9 X 10^ GeV^ (^) = ^ ^^^^^ 



(44) 



with C given by Eq. ()42|1 . The neutrino mass is quadratic in C, and thus the dependence 
on the uncertain parameters grouped in C is larger than in case 1. Minimizing C gives 
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Mi > 2 X 10^3 GeV. 

The lower bound on the neutrino mass Eq. is shown in Fig. Q as well as the 
kinematic constraint vn^j2 < < ~ r7/2, and the gravitino constraint. The reheat 

temperature is independent of the neutrino mass, but does depend on the coupling k. For 
K < 10^^ one has Tr < 3 x 10^ GeV and there is no gravitino problem. Leptogenesis is 
possible for Mj ~ 10^^ — 10^^ GeV and k ~ 10^^. The RH neutrino mass increases rapidly 
with small and much smaller couplings are excluded. The large neutrino masses needed 
are incompatible with a non-renormalisable mass term as in Eq. (fTT|) and perturbative 
couplings, see Eq. (j20|l . 

4.3.3 Case 3 

Consider now the case that the lightest RH neutrino reaches thermal equilibrium after 
inflation (Mi < Tr). As follows from Eq. (|TIHl thermal equilibrium can only occur if the 
infiaton decays in N2 or N^, and not in the lightest RH neutrino. The lepton asymmetry 
is thus produced via string decay in one of the heavier RH neutrinos and i = 2, 3. Any 
produced lepton asymmetry will be washed out until L-violating reactions fall out of 
equilibrium at T ~ Mi. We assume that this occurs in the scaling regime. 

Loop scenario The lepton number density is now given by Eq. (p?T|l . Diving by the 

entropy we get 



s 'jla \mp 



' in 



3/2 



(45) 



with ifin ~ Mf/rrip. Setting it equal to the observed value gives the RH neutrino mass 
needed for successful leptogenesis: 

M, ^ 10^^ GeVC'/' l^j > 2 X 10^^ GeV ( ^^^^ j . (46) 

In the second step we used the equilibrium condition Eq. (fT^ and r] = 77^^; further we 
defined 

C = (47) 

V In 

The number of RH neutrinos released per loop is bounded by xa? < (see Eq. fjHIH) ). 
Further a gives the loop size at birth (see Eq. P^ ): taking a = (FG/x)" with n = 3/2 
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Figure 5: Mj vs. k for case 3 with rj = rj^^, Af = I, and a = ai. The hla part corresponds 
to leptogenesis in the loop scenario with 1 < a;jv < and the hnes A-B to the VHS 
scenario with A = {•ji^, X , f^^^) = (7s, 0,1), and B = {%,A,1). The parameter space is 
bounded by the kinematic constraint (top), and the gravitino constraint (parallel lines, 
forTR = 10i\...,10i^GeV). 

instead of n = 1 lowers C by about a factor 10. Thus C > O.lfi;. The lower bound on 
the RH neutrino mass is then 



M, > 6 X 10^2 GeVAT^/^ ( 



K 



10- 



1/4 



(48) 



K 



-2 



together with the CMB constraint 10"^ < ^ < lO^V-^- 

The bound on the neutrino mass Eq. pUj) is shown in Fig. (jSJ for 1 < /^v < 
together with the kinetic and gravitino constraint. Large couplings k > 10~^ are needed, 
which is marginally excluded by the CMB data. The reheat temperature has to be large 
T^> lOi^GeV. 



VHS Scenario The lepton number density is now given by Eq. (jS^), leading to 



nL_ 

s 



3/2 

7in mxmp 
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(49) 



with i/in ~ Mf/nip, = 7^ ~ 1, and with i = 2,3. Leptogenesis requires 



.m, GeV^/2 / M 7 x lO^^ GeVC 



where in the second step we used the equihbrium condition Eq. (fTIJ|) and r] = r]^^. 

The EGRET data requires the fraction of the total energy that goes into X-particles to 
be small < 10~^; as a result leptogenesis is not efficient enough and too large neutrino 
masses are needed, incompatible with the kinematic constraint. If however in the early 
scaling regime f^^l also possible — remember that the EGRET bound is determined 
by late times — smaller neutrino masses are possible, as shown in Fig. (jSJ. 

Our results Eq. (pHt lKU jl for C = 1, C" = 1 agree with those found in Ref. (THj. Ref. [Ej 
assumes degenerate light neutrino masses, and finds stronger bounds. 



5 /i-term 

The /i-problem can naturally be resolved in SUSY hybrid infiation with the introduction 
of a superpotential term 

W = XSHH' (51) 

where H, H' contains the two Higgs doublets of the MSSM. After infiation, S gets a 
VEV due to low energy SUSY breaking which is of order {S) ~ ?Ti3/2//« provided A > /t 
(otherwise 5* ends up in wrong minimum), and the /i-term is generated. However, the 
superpotential term above opens up a new decay channel for the infiaton and jeopardizes 
non-thermal leptogenesis. 

The kinematic constraint Mj < together with the constraint \ > k assures that 

decay rate for infiaton decay into SM higgses and higgsinos 

A2 



H 



-m^ (52) 



167r 

is larger than the decay rate into RH neutrinos. Infiaton decay is predominantly into SM 
Higgs fields (unless all the couplings are tuned 2Mj/?7 ~ k ~ A), and NT leptogenesis via 
infiaton decay does not occur. Likewise, string decay into Higgs fields does not contribute 
to the lepton asymmetry, since the Higgs decays into SM Higgses. On the other hand, 
if the string decays into gauge fields, or if RH zero modes are released during decay, 
leptogenesis is still possible. We will consider this possibility in some detail. 
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As a side remark, we note that there are other ways in these models to generate a 
/i-term which do not alter the inflaton decay rate, and are thus compatible with NT 
leptogenesis from reheating j38j . 

Either case 2 or case 3 is realized, depending on whether the reheat temperature 



is smaller or larger than the mass of the lightest RH neutrino. Note that the reheat 
temperature is minimized in the limit X ^ k. 

5.1 Results 
5.1.1 Case 2 

Leptogenesis is a result of string decay into RH neutrinos and into gauge fields which in 
turn decay into RH neutrinos. The RH neutrinos are out-of-equilibrium at all times, and 
the asymmetry is dominated by the initial time just at the end of inflation. There is no 
contribution to the asymmetry from inflaton decay nor from thermal leptogenesis. 

The lepton asymmetry is given by Eq. (jHS)) with the replacement Tn ^ Th and 
X = A. This gives 



where we used rj = rj^^. The reheat temperature is a function of k only: Tr(k/A) = 
10^ 10^°, 10^2 GeV for k = 10-^ 2 x 10"^, 3 x lO^^. The neutrino mass can be lowered 
by lowering the ratio (k/A) < 1, but at the cost of increasing the reheat temperature 
with the inverse ratio. This makes it harder to satisfy the out-of-equilibrium condition 



The results are shown in Fig. IHl for various values of C together with the kinematic 
constraint Mj < ~ t], and the out-of-equilibrium condition Mj > Tr. Note that the 
lower bound on the neutrino mass is proportional to f~^, and the results for different 
values of can be obtained by multiplying with the appropriate factor. Masses as low 
as Mi ~ 10^ GeV are compatible with leptogenesis. 




(53) 




(54) 



Mi > Tr. 



5.1.2 Case 3 



The lightest RH neutrino reaches thermal equilibrium, and all asymmetry is erased until 
it falls out of equilibrium at T ~ Mi. 
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Figure 6: Mj vs. n for case 2 with /x-term and rj = M = 1. The hnes A-B correspond 
to A = (7in, X, f^) = (7min5 ^5 1), and B = (7best5 A, !)• The parameter space is bounded 
by the kinematic constraint (top), and the wash-out constraint Tr < Mi (bottom). 

Loop scenario The lepton asymmetry is independent of the decay rate, in particular, 
on whether the inflaton decays into RH neutrinos or SM Higgses. Hence, the results of 
case 3 without a /i-term apply and 



The only differences are that now z = 1 is possible and still the lightest RH neutrino 



decay rate and the symmetry breaking scale rj; the K-dependence enters only via xn- 

The lower bound on Mi is shown in Fig. [3 for 1 < xn < together with the 
kinetic and equilibrium Mi < Tr constraint. Here it is assumed that (/t/A) = 1, which 
minimizes the reheat temperature. Hence, the equilibrium constraint can be relaxed by 
taking (k/A) < 1, but at the cost of increasing the reheat temperature, and thereby 
aggravating the gravitino problem. 




(55) 



reaches equilibrium, and /at < k ^ so that C > k 



The results do not depend on the 
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Figure 7: Mj vs. /t for case 1 with /i-term and rj = rj^^, A/" = 1. The hla part corresponds 
to leptogenesis in the standard scenario with 1 < /^v < k^^, and the hne A to the VHS 
scenario with (7in, X, f^) = 0, 1). The parameter space is bounded by the kinematic 
constraint (top), and the equihbrium constraint Tr > Mi (top). 



VHS scenario The lepton asymmetry is the same as in case 3 without a /i-term, and 
requires 

M ~ 4 X lO'Uhi^ (Mi) - ^xlO-^GeVC- ( M,\ 

M._4xiu y^^j - i^jj^^^ y^j (50) 

where in the second step we used t] = rj^^. The result is shown in Fig. [7| The lower 
bound on the neutrino mass is proportional C oc /~^^^, and the result for different values 
of can be obtained by multiplying with the appropriate factor. 



For couplings k, < 10~^ the reheat temperature is Tr < 10^° GeV, and this is a scenario 
with can accommodate inflation and leptogenesis, and in which both the gravitino problem 
and the /i-problem are solved. The lightest RH neutrino can be in equilibrium or not: 
both case 2 and 3 can work. 
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6 Conclusions 



In this paper, we have investigated various possibihties for leptogenesis after hybrid infla- 
tion when gauged B-L is spontaneously broken at the end. One of the Higgs fields gives 
heavy Majorana mass to the RH neutrinos and NT leptogenesis can take place during 
reheating via infiaton decay into RH (s)neutrinos. Cosmic strings form at the end of 
inflation [HI EI- If stable, they also contribute to primordial fluctuations. Interestingly 
enough, since the string Higgs field breaks these are the so-called B-L strings whose 
decay gives a second NT contribution to the lepton asymmetry of the universe In 
this paper we investigated which of these two mechanisms is most efficient, taking into 
account the CMB constraints [T3] . 

Leptogenesis via infiaton decay can account for the observed asymmetry for neutrino 
masses in the range Mj = 10^ — 10^^ GeV, and quartic Higgs couplings n = 10^^ — 10~^. 
The minimal reheat temperature required for getting enough lepton asymmetry is Tr ~ 
7 X 10^ GeV. To assure that all are out-of-equilibrium at Tr, and there is no wash-out 
of lepton number, the mass of the RH neutrino the infiaton decays into has to satisfy 
M, < 102 X Tn. 

The calculation of the lepton asymmetry created by string decay is hampered by our 
poor knowledge of the properties of the string network. The initial string density, the loop 
formation and decay mechanisms all introduce uncertainties. In general, only the 'best 
case' scenarios give a large contribution to the lepton asymmetry, in which it is assumed 
that the initial string density is high and/or that the fraction of string energy going into 
X-particles is appreciable. We argued that both of these assumptions are not far-fetched, 
as they fit well with our knowledge of cosmic strings. 

We distinguished three different cases, depending on whether the infiaton field decays 
into RH (s)neutrinos, and wether the lightest neutrino is out-of-equilibrium at reheating. 

In case 1, the infiaton decays into RH neutrinos and all RH neutrinos are out-of- 
equilibrium at Tr. The lepton asymmetry is determined by the energy density in the 
string network right at the end of inflation. If the strings mostly decay into Higgs particles, 
the contribution to the asymmetry is subdominant with respect to the contribution from 
infiaton decay. The strings can give a dominant contribution if it decays mostly into RH 
neutrinos or into gauge particles, but for the latter only in the 'best case' scenario. 

If infiaton decay is not into RH neutrinos and the RH neutrinos never attain ther- 
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mal equilibrium (case 2), the only contribution to the baryon asymmetry of the Uni- 
verse is from cosmic strings decay. String decay into RH neutrinos and vector particles, 
which subsequently decay into RH neutrinos, can still produce a lepton asymmetry. The 
minimal reheat temperature required for getting enough lepton asymmetry Tr ~ 10^ 
GeV. Gravitational infiaton decay can ameliorate the gravitino constraint. However, large 
Mi ~ 10^^ — 10^^ GeV RH neutrino masses and couplings k > 10"'^ are required. Such 
large couplings favor small Higgs representations. 

The most unfavorable scenario is case 3, in which the lightest RH neutrino attains 
thermal equilibrium. The mass of the heaviest RH neutrino the infiaton can decay into 
(which cannot be Mi) must be large. Mi > 10^^ GcV, the coupling must also be large, 
and the reheat temperature must be high Tr > 10^'^ GcV. We note that in this case there 
is also a contribution to the lepton asymmetry from thermal leptogenesis. 

Finally, we looked at the possibility of generating the MSSM /x-term dynamically. 
Infiaton decay is now into SM Higgses and Higgsinos, and only case 2 and 3 can occur. 
For a wide range of RH neutrino masses and for couplings k < 10~^ this is a scenario which 
can accommodate inflation and leptogenesis, and in which both the gravitino problem and 
the //-problem are solved. 
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A CP- asymmetry 



In this appendix we derive the typical values for the CP-asymmetry factor e^, extending 
the usual analysis to i = 2, 3. We use the formalism and notation of [5^ . 
The CP-asymmetry is 

1 1 

^.=-i^pt)7i:'™[{(«'%n/(ii4) (57) 

with in SUSY theories 

2 



X 



+ ln(l + 



(58) 



For hierarchical RH neutrino masses Mi -C M2 <^ M3, we need the limits a; <^ 1 (where 
/ —>■ 3/^/x) and a; ^ 1 (where / —^fx{2 + ln(2;))). The function / is maximized for 
f ^ 1, but this does not occur for hierarchical RH masses. Then 

3 1 



ei 



11 



Im[{(AAt),,n^ + Im[{(AAt),3p]^ 



' <I. (59) 



Here / ~ 0{1) a phase factor, v = 174GeVsin/3 the Higgs VEV (note that sin/3 pa 1 for 
/? > 3). In the second step we have used 

{hh%2 ~ {ml/v^)l22, {hh^)i2 ~ {mumjv^)li2, {hh%, ^ {mtmjv^)l23, (60) 

with I order one constants, and mu,mc,mt he Dirac neutrino masses which are labeled 
in analogy with the quark masses. One can express the Dirac masses in terms of the RH 
neutrino masses: 

\Mi\=ml/Ai, \M2\=ml/A2, \Ms\ = m^JA^, (61) 



with Ai = 5^2^^"^ Li' ^2 = \/^"^ltm/2 and A3 = 2|mi|/s^2- The numerical values for 
Ai depend on the spectrum of light neutrino masses. For example Ai = 2x 10^^^ GeV, 3 x 
10~^^ GeV for normal hierarchy (ma ^ (Am'l^^Y^'^ ^ m2 ~ (Am^Qj)^/^ ^ mi) respec- 
tively inverted hierarchy(mi ^ m2 ~ (Am^^^^)^/^ ^ ^ (Am^^j)^/^) [39j. Here it is 
approximated that ^^13 = and 6'23 = 7r/4. The CP-asymmetry is 

3AiMi ^.,f Kh 



V 106 GeV.' ^ ' 
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The bound in Eq. ()22|) is obtained for inverted hierarchy and order one phases. For quasi- 
degenerate neutrinos (mi ^ m2 ~ ms ^ (Am^^jj^)^/^) the asymmetry is suppressed by a 
factor Ay^mf^/m with rh = + ml + m| fT7\ . 

A similar calculation for decay of the next to lightest RH neutrino gives 
1 1 



£2 



in {hW) 
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Im [{(/./.%}2]^ + Im [{(/./.t),3f ] ^ 



1 

2 X 10 



ml Ml A3 



-10 



M2 
M2 



(63) 



,106 GeVy 

where in the last step we have neglected the subdominant term proportional to Mi, and 
we have used ^3 < 6 x 10^^^ GeV valid for hierarchical (light) neutrinos. The bound is 
of the same order of magnitude as the bound on ei, but with Mi ^ M2. Finally, for the 
CP-asymmetry of the heaviest RH neutrino we get 



MfAi M|A2 



+ 



3-7x 10 



-11 



Ms 
M2 



A2M2 M2 
Sttv^ M3 



M2 
Ms 



(64) 



^106 GeV, 

5 X 10^^^ valid for hierarchical light neutrinos. Note that 



where we have used A2 = 3 
es is suppressed by a factor M2/M3 and is smaller than ei, 62 
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